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Condensation phase transitions of symmetric conserved-mass aggregation model on 

complex networks 
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We investigate condensation phase transitions of symmetric conserved-mass aggregation (SCA) 
model on random networks (RNs) and scale-free networks (SFNs) with degree distribution P(k) ~ 
fc -7 . In SCA model, masses diffuse with unite rate, and unit mass chips off from mass with rate ui. 
The dynamics conserves total mass density p. In the steady state, on RNs and SFNs with 7 > 3 for 
uj 7^ 00, we numerically show that SCA model undergoes the same type condensation transitions 
as those on regular lattices. However the critical line pc(u>) depends on network structures. On 
SFNs with 7 < 3, the fluid phase of exponential mass distribution completely disappears and no 
phase transitions occurs. Instead, the condensation with exponentially decaying background mass 
distribution always takes place for any non-zero density. For the existence of the condensed phase 
for 7 < 3 at the zero density limit, we investigate one lamb-lion problem on RNs and SFNs. We 
numerically show that a lamb survives indefinitely with finite survival probability on RNs and SFNs 
with 7 > 3, and dies out exponentially on SFNs with 7 < 3. The finite life time of a lamb on SFNs 
with 7 < 3 ensures the existence of the condensation at the zero density limit on SFNs with 7 < 3 
at which direct numerical simulations are practically impossible. At ui — 00, we numerically confirm 
that complete condensation takes place for any p > on RNs. Together with the recent study on 
SFNs, the complete condensation always occurs on both RNs and SFNs in zero range process with 
constant hopping rate. 

PACS numbers: 05.70.Fh,05.40.-a,89.75.Da,89.75.Hc 



I. INTRODUCTION 

Nonequilibrium condensation phase transitions from 
fluid phase into condensed phase have been observed in 
a variety of phenomena rangingfrom traffic flow to poly- 
mer gels 0, 1, S 1 1 1 Hi, I E [H E El. In the 
steady state, a finite fraction of total particles condenses 
on a single site in condensed phase when the total parti- 
cle density p is increased beyond a certain critical value 
p c . In fluid phase below p Cl the particle number of each 
site fluctuates around p without the condensation. 

Various kinds of nonequilibrium mass transport mod- 
els exhibit the condensation transitions or only conden- 
sation. The simplest and well-known model is zero-range 
process(ZRP) in one-dimension p). In ZRP, many iden- 
tical particles occupy sites on a lattice. Each site may 
contain an integer number of particles and one of these 
particles can hop to one of the nearest neighboring sites 
with a rate that depends on the number of particles at the 
site of departure. The chipping (single-particle dissocia- 
tion) and aggregation processes of ZRP describe various 
condensations such as jamming of traffic 0, bunching of 
buses , coalescence of shaken steel balls Q and conden- 
sation of edges in networks : 5j. Recent studies of ZRP 
on scale-free and directed networks reveal the con- 
ditions under which condensation takes place. 

Another important class of condensation transitions 
emerges when the diffusion of the whole particles of a 



single site is involved in addition to the chipping and ag- 
gregation. These processes arise in a variety of phenom- 
ena such as polymer gels @, the formation of colloidal 
suspensions jl| , river networks 0, E] and clouds E| . 
Recently studied conserved-mass aggregation(CA) model 
is the simplest one incorporating^ diffusion, chipping and 
aggregation upon contact 0, I n one dimensional 

CA model, the mass m, of site i moves either to site 
i — 1 or to site i + 1 with unit rate, and then m, — > and 
TOi-ti — > my-i+mj. With rate to, unit mass chips off from 
site i and moves to one of the nearest neighboring sites; 
m j — > m,i — 1 and mm — > mj±i + 1. As total masses are 
conserved, the conserved density p and u> determine the 
phase of CA model. The condensation transition arises 
via the competition between diffusion and chipping pro- 
cesses. The diffusion of masses tends to produce massive 
aggregates and consequently creates more vacant sites. 
The chipping of unit mass tends to prevent the forma- 
tion of aggregates, so that it leads to a replenishment of 
the lower end of mass distribution. 

The single site mass distribution P(m), i.e., the prob- 
ability that a site has mass to in the steady state, was 
shown to undergo phase transitions on regular lattices 
[l3| . For a fixed lu, as p is varied across the critical value 
jfc^ Q; f ne behavior of P(m) for large to was found to be 

- m/m " p < PcH, 

P{m)~{ m - P = Pc(u), (1) 

<r T + infinite aggregate p > p c (uj), 
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where p c is given as p c (w) = \/uT+T — 1. p c and r are 
shown to be independent of spatial dimension d [TJ| . The 
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tail of the mass distribution changes from exponential to 
an algebraic decay as p approaches p c from below. As 
one further increases p beyond p c , this asymptotic alge- 
braic part of the critical distribution remains unchanged 
but in addition an infinite aggregate forms. This means 
that all the additional mass (p — p c )L d condenses onto a 
single site and does not disturb the background critical 
distribution. The u> — oo case corresponds to ZRP with 
constant chipping rate, and then there is no condensation 
transitions on regular lattices. The critical exponent t is 
same everywhere on the critical line p c {ou). Recent stud- 
ies showed that r depends on the symmetry of movement 
and constraints of diffusion rate [jjj, . In what follow- 
ing, we only consider the symmetric CA (SCA) model 
where diffusion and chipping directions are unbiased |l3j | . 

As mean field theory correctly predicts the phase di- 
agram and r of SCA model in one dimension 0], the 
nature of the condensation transitions should remain un- 
changed in any higher dimensions. However recent stud- 
ies of various dynamics on network structures showed 
that the structure of networks leads to more rich and 
intriguing behavior different from that predicted by the 
standard mean-field theory on regular lattices [a, IfJ 113 • 
For example, ZRP on SFNs was shown to exhibit conden- 
sation even at constant hopping rate which is the case of 
an infinite chipping rate of SCA model [|| . 

In this paper, we investigate the effect of network struc- 
tures on the condensation transitions of SCA model using 
random networks (RNs) and scale-free networks (SFNs). 
As we shall see, on RNs and SFNs with the degree expo- 
nent 7 > 3, SCA model undergoes the same type of con- 
densation transitions as those in regular lattice across a 
critical line p c (w) in p-uj plane with the exponent r = 5/2. 
However on SFNs with 7 < 3 where one or several nodes, 
so-called hub nodes, have a finite fraction of links, the 
fluid phase completely disappears and the condensation 
with exponentially decaying background mass distribu- 
tion takes place for any non-zero density. The out-line 
of this paper is as follows. In Sec. II, we introduce SCA 
model on complex networks. The condensation transi- 
tions on RNs and SFNs are discussed in Sec. Ill and IV. 
To understand the condensation on SFNs with 7 < 3, we 
discuss lamb-lion problems on SFNs in Sec. V. In Sec. 
VI and VII, we discuss SCA model at ui — 00 and the 
effect of diffusion of masses on average mass distribution 
on degrees respectively. Finally we summarize our results 
in Sec. VIII. 



II. SCA MODEL ON NETWORKS 

We consider a network with N nodes and K links. 
The degree ki of a node i is defined as the number of its 
links connected to other nodes. The average degree of 
a node < k > is given as < k >= 2K/N. The degree 
distribution P(k) is a Poisson distribution for RNs and 
a power-law distribution of P{k) ~ fc~ 7 for SFNs. Each 
node may have an integer number of particles, and the 



mass of a node is defined as the number of particles at 
the node. Initially M particles randomly distribute on 
N nodes with given conserved density p = M/N . Next 
a node i is chosen at random and one of the following 
events is occurred: 

(i) Diffusion: With unit rate, the mass rrii moves to 
the randomly selected nearest neighboring node j . If the 
node j has already mass rrij , then the aggregation takes 
place; m, — > and rrij — > rrij + mi. 

(ii) Chipping: With rate oj, unit mass (a single parti- 
cle) at node i chips off and moves to the randomly se- 
lected nearest neighboring node j; m, — * rrii — 1 and 
rrij —>■ rrij + 1. 

As mentioned in Sec. I, diffusion and chipping pro- 
cesses compete each other so a condensation transition 
physically depends on two external parameters p and lj. 
At u> = 0, however, only diffusion and aggregation take 
place so masses always condense onto a single site (com- 
plete condensation). On the other hand, the lo = 00 case 
corresponds to ZRP with constant chipping rate. It was 
shown that ZRP on scale-free networks exhibits complete 
condensation when mass-dependent chipping rate u(m) 
is given as u(m) ~ m 5 , but <5 < 6 C , where S c is some 
threshold value [f|. Hence the SCA model of u> = 00 is 
just the S — case of ZRP so the complete condensation 
always aries at lo = 00. As we shall see in the following 
section, the complete condensation also takes place on 
random networks at p > 0. 

For the construction of SFN, we use a static model |l8| 
instead of preferential attachment algorithm |5j] . In static 
model, it is desired to use large < k > to construct fully 
connected networks. In simulations, we use < k >= 4. 

III. SCA MODEL ON RNS AND SFNS WITH 

7 > 3 

We consider SCA model on RNs and SFNs with 7 > 3 
having nodes N = 10 4 and links K — 2 x 10 4 for w^oo. 
We define the single node mass distribution P(m, t) at 
time t, i.e., the probability that a node has mass m at t 
as 

1 1 N 

Pim, t) = — n(m, t) = — ^ <f mi (t),m . ( 2 ) 

i=l 

where n(m, t) is the number of nodes with mass m at 
time t and rrii{i) is the mass of node i at t. From the 
normalization condition of P(m, t), i.e., Yl m P{ m i t) = 1, 
all P(m, t) should reach the stationary state at the same 
time. For instance, P(m = 0, t) is given as P(0, t) = 
1— J2 m ^o P( m > Hence, if P(0.t) saturates to its steady 
state value P(0) after some time tq, then P{m, t) for any 
mass m also should reach its steady state value P(m) 
after the time To . So the characteristic time r m at which 
P(m, t) saturates is expected to be independent of mass 
m, and 

to = n = ••• = T m * , (3) 
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where m* is the maximal mass. In Fig. 1, we plot P(m, t) 
of m = 1 and 10 for u> = 1.0 and p = 0.1 on various net- 
works. The vertical dotted line denotes the characteristic 
time T m , which supports Eq. 

We measure the steady state distribution P(m) by av- 
eraging P(m, t) after r m . Fig. 2 shows the plots of P(m) 
versus m for RN and SFN with 7 = 4.3 with u> = 1.0. 
The critical density p c is p c (oj = 1) = 0.31(3) for RN and 
0.24(3) for SFN with 7 = 4.3 respectively. In fluid phase 
(p < p c ), P(m) exponentially decays for large mass in 
both RN and SFN with 7 = 4.3. At p = p c , P(m) alge- 
braically decays on both RN and SFN with 7 = 4.3 as 
P{m) ~ m~ T with t = 2.38(3) for RN and r = 2.33(2) 
for SFN with 7 = 4.3 without condensations. The values 
of t of both networks also agree well with r = 5/2 of reg- 
ular lattices 00]. In condensed phase (p = 3.0 > p c ), 
the access mass of (p — p c )N « 2.7 x 10 4 condenses on 
a single node in the both networks without changing the 
background critical distribution. To map out phase dia- 
grams for RN and SFN with 7 = 4.3, we measure p c and 
t for oj = 0.1, 1.0 and 10, and find that t is the same 
within error on the critical line. Fig. 3 shows the phase 
diagrams of RN and SFN with 7 = 4.3. The solid line 
denotes the critical line, p c — y/uT+T — 1 of regular lat- 
tices 0. On regular lattices, the p c (oj) is independent of 
dimensionality. However on networks, p c {^>) depends on 
the underlying network structures. We also confirm that 
SFN with 7 = 3.5 undergoes the same type of condensa- 
tion transitions. Therefore we conclude that SCA model 
on RNs and SFNs with 7 > 3 undergoes the same type of 
condensation transitions from fluid phase into condensed 
phase for < lu < 00 as those on one-dimensional regular 
lattice. However the critical line depends on the under- 
lying network structures. 



coexist on SFNs with 7 < 3. Such a behavior was found 
in the lattice gas model proposed to describe the distri- 
bution of droplets in the fragmentation process following 
a nuclear collision • 

The nonexistence of the fluid phase for 7 < 3 mainly 
comes from the hub structure of SFNs with 7 < 3 where 
one or several nodes have a finite fraction of total links 
|5j. Due to the hub structure, all masses move to a 
hub node by both diffusion and chipping processes un- 
like SFNs with 7 > 3 where the chipping tends to split 
masses and prevents aggregations. Hence there are no 
processes to prevent the formation of infinite aggregation 
for 7 < 3. The difference from ZRP is that the diffusion 
process moves an infinite aggregation of mass rrihub at a 
hub node to the others. Then the chipping distributes 
small masses of rrihub onto neighboring nodes. In this 
way, nodes with small degree can have mass so the resul- 
tant P(m) is exponential. 

It is hard to show the existence of infinite aggregation 
for the limit of p — ► numerically with finite size net- 
works. Instead we investigate annihilating random walks 
of two particles on SFNs (lamb- lion problem) for the con- 
dition of the existence of the fluid phase on SFNs. For 
p — > limit, we assume an infinite aggregation and ne- 
glect the presence of other masses. Then unit mass chips 
off from the aggregation, and diffuses around on the net- 
works. For the existence of an infinite condensate in the 
steady state, the two masses should aggregate again in 
finite time interval. If not, unit mass continuously chips 
off from the infinite aggregation, which will finally disap- 
pear. Hence if the unit mass (lamb) survives indefinitely 
without meeting the infinite aggregation (lion), then the 
fluid phase exists at zero density limit. Otherwise, the 
condensed phase exists for any p > 0. In next section, 
we investigate a lamb-lion problem on RNs and SFNs. 



IV. SCA MODEL ON SFN WITH 2 < 7 < 3 



LAMB-LION PROBLEM ON NETWORKS 



For oj 7^ 00, we perform simulations on SFNs with 
7 = 2.4 and 3 with N = 10 4 nodes and K = 2 x 10 4 links. 
We plot P(m, t) of 7 = 2.4 in Fig. 1 (c) for example, 
and average P(m,t) after r m for the steady state P(m). 
As mentioned in the previous section, r m is same for all 
masses. 

The mass distribution P(m) for 7 < 3 show quite dif- 
ferent behavior from that for 7 > 3. Fig. 4 shows P(m) 
on SFNs with 7 = 2.4 and 3 for p = 0.2 and 3.0 with 
u> = 1. As in the condensed phase of RNs and SFNs 
with 7 > 3, there is a condensation of mass m* pN . 
However P(m) exponentially decays for large mass rather 
than power-law. We also measure P(m) for various to 
and p, and confirm the same behavior of P(m). There- 
fore we conclude that for w ^ 00, SCA on SFNs with 
7 < 3 exhibits an infinite aggregation with an exponen- 
tial background mass distribution. We call this phase 
incomplete condensed phase to reflect exponential mass 
distributions (Fig. 4 (c)). Intriguingly the two features 
of the fluid and the condensed phase on regular lattice 



We consider one lamb and one lion problem on net- 
works. If a lamb meets a lion on the same node, then it 
dies. Hence the interesting quantity is the survival prob- 
ability S(t) of a lamb at time t. On d dimensional regular 
lattice, S{t) depends on d as follows |23ll24l l25ll26| . 



' t d/2-l fo r rf<2 

S(t)~\l/lnt ford = 2 
constant for d > 2 



(4) 



On networks, the steady state distribution of a random 
walker is proportional to the degree ki of node i 0] as 



N 



(•5) 



Hence the probability of finding two walkers at the same 
node should depend on the second moment < k 2 >. 
While < k 2 > is finite for RNs and SFNs with 7 > 3, 
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it diverges for SFNs with 7 < 3. So RNs and SFNs 
with 7 > 3 correspond to infinite dimensional homoge- 
neous regular lattice in macroscopic scale, while SFNs 
with 7 < 3 are generically inhomogeneous. Therefore 
S(t) is expected to be constant in the thermodynamic 
limit on RNs and SFNs with 7 > 3. It means the ex- 
istence of the fluid phase for the limit of p — for any 
lo 7^ 00 as shown in Fig. 3. On the other hand, for the 
incomplete condensed phase of 7 < 3, S(t) should decay 
to zero for 7 < 3. 

We measure S(t) of a lamb on RNs and SFNs with 
various 7 for several values of N from 10 3 to 10 6 and up 
to 10 6 time steps. Initially a lamb and a lion are placed 
on randomly selected two nodes. In Fig. 5, we plot S(t) 
for 7 = 4.3 (a) and 2.4 (b) for various system size N. 
Interestingly S(t) decays exponentially in both SFNs for 
any N. For example, S(t) of 7 = 4.3 and N = 10 6 ex- 
ponentially decays with very large characteristic time r 
order of 10 6 over the whole time interval. We also con- 
firm the same exponential decay of S(t) on RNs for any 
N. Unlike on regular lattices where the exponential de- 
cay of S(t) comes from finite size effect of lattices, S(t) 
on networks shows such a exponential decay at very early 
time. The exponential decay of S(t) mainly comes from 
the small world nature of networks i.e., the average 
path length between a pair of nodes increases logarith- 
mically with N. As S(t) shows the exponential decay for 
any 7, we are interested in the average life time T of a 
lamb rather than S(t) itself. 

We measure the life time T of a lamb on RNs and SFNs 
with several 7 values for N up to 10 6 and plot T against 
N in Fig. 5 (c) and (d). In Fig. 5 (c), T linearly increases 
with N for both RN and SFN with 7 = 4.3. Assuming 
T ~ N a : we estimate a = 1.0(2) for RN and SFN with 
7 = 4.3 by measuring slopes. However for 7 < 3 of Fig. 
5 (d), T exhibits quit different behavior. T of 7 = 2.15 
and 2.4 tend to saturate to some asymptotic value with 
decreasing successive slopes. However the characteristic 
size N s ("f) at which T begins to saturate increases as 7 
approaches 7 = 3. For example, T of 7 = 2.75 seems to 
algebraically increase up to N = 10 6 . It means that N s of 
7 = 2.75 is already larger than 10 6 . To see the saturation 
at 7 = 2.75, big size N much larger than 10 6 is needed but 
it is practically very difficult. It is hard to believe that T 
diverges in power-law fashion as 7 approaches 3 because 
T of small 7 is already begin to saturate at moderate N 
order of 10 5 ~ 10 6 . Hence T of 7 < 3 is expected to 
saturate to its asymptotic value even though we cannot 
see the saturation via simulations as 7 approaches 3 from 
below. At 7 = 3, T also increases algebraically with N. 
By measuring successive slopes, we estimate a = 0.85(1) 
for 7 = 3. However, as the exponent a continuously 
varies as 7 increases to 3, the value of a at 7 = 3 may 
have no special meaning. As shown in Fig. 4 (a) for 
7 = 3, the condensation already occurs at sufficiently 
low density p — 0.2 even for small size N = 10 4 so T 
of 7 = 3 is believed to saturate in thermodynamic limit 
as for 7 = 2.75. Therefore we are convinced that the 



life time of a lamb is finite in thermodynamic limit of 
N — > 00 for 7 < 3. From the behavior of the life time 
T, the asymptotic behavior of S(t) is investigated in the 
following. 

As shown in Fig. 5 (a), S(t) decays exponentially 
as S(t) = S e~ t/T . From the definition of T = 
J °° t(—dS(t)/dt)dt, we have T ~ r. Hence r diverges 
linearly on SFNs with 7 > 3 and RNs, while it saturates 
for 7 < 3. The scaling behavior of r means that S(t) on 
RNs and SFNs with 7 > 3 is nonzero constant in ther- 
modynamic limit of N — > 00, and exponentially decays 
to zero on SFNs with 7 < 3. Hence we have 

BmmJ^ ;^ 3 j (6) 
v ; ISoo ( 7 >3) , K ' 

where Too is the asymptotic value of r, and depends on 
7. As S(t) decays exponentially to zero for 7 < 3, two 
particles initially located on the same node should meet 
or aggregate within finite time interval. Hence the finite 
life time on SFNs with 7 < 3 supports the numerically 
expected condensed phase for any p > on SFNs with 
2 < 7 < 3 as mentioned in the previous section. On the 
other hand, the finite Soo on RNs and SFNs with 7 > 3 
implies the existence of the fluid phase at p — > limit, 
which is also consistent with the phase diagrams of Fig. 
3. 



VI. SCA MODEL AT u = 00 

At lo = 00, SCA model corresponds to ZRP with con- 
stant hopping rate on networks. On SFNs with 7 > 2, it 
was shown that ZRP undergoes complete condensation 
for any p > [|| . In condensation phenomena of ZRP on 
SFNs, the average mass at a node with degree k, in 
the steady state exhibits a discontinuous jump at a hub 
node. In ZRP with constant chipping rate for example, 
mfc linearly increase until k < khub and jumps to the 
value m* w pN at k = khub |6|- The khub is the degree 
of the hub node which has the largest number of links. 

To see whether the condensation occurs on RNs in 
ZRP, we measure rrik for p = 0.4 and 3.0 on RNs with 
N = 10 4 . As shown in the bottom inset of Fig. 6 (a), 
rrifc increases linearly in k and finally jumps at khub as 
in ZRP on SFNs. For a fixed N, the mass of khub, rrihub 
decreases as p goes to zero as shown in the inset of Fig. 
6 (a). So it is hard to observe the condensation at low 
density for small TV. Since the complete condensation of 
ZRP on networks mainly comes from the inhomogeneity 
of degree distribution rather than the competition be- 
tween dynamical processes, ZRP on RNs should exhibit 
the complete condensation for any p > as on SFNs. At 
lo = 00, we also measure on SFNs with 7 < 3 and 
confirm the complete condensation as shown in Ref. 0. 
Together with the recent study on SFNs 6], the com- 
plete condensation always occurs on both RNs and SFNs 
in zero range process with constant hopping rate. 
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The discontinuity of rrife clearly comes from the na- 
ture of chipping process of unit mass. That is, unit mass 
hops around according to the steady state distribution 
of Eq. (5) , and then there is enough net particle current 
into the hub nodes to form infinite aggregation. However 
when the diffusion of the whole mass of each node turns 
on, the aggregation at the hub node diffuses around all 
nodes. The diffusion of the condensate completely change 
the behavior of mk in the condensed phase of finite lo. mk 
linearly increase up to khub without any jump in the con- 
densed phase of lo ^ oo. In next section, we numerically 
and analytically confirm the linearity of mk in k. 



VII. AVERAGE MASS OF A NODE WITH 
DEGREE K 



Another interesting quantity in condensation phenom- 
ena on networks is the average mass at a node with de- 
gree k, rrifc in the steady state. In ZRP with chipping 
rate u(m) ~ m 5 on SFNs Q, the complete condensation 
takes place for S < S c = 1/(7 — 1). For S < 8 C , rat linearly 
increases with k for k << k Cl and algebraically increases 
as fc 1 / 5 for k > k c . Especially at 6 = 0, m k linearly in- 
creases until k < khub and jumps to the value m* w pN . 
The khub is the degree of the hub node which has the 
largest number of links. The same type of behavior on 
RN for lo = 00 is shown in the inset of Fig. 6 (a). 

In SCA model, as all masses can perform random walks 
according to the steady state distribution P°° of Eq. (5), 
one can expect the jump of m k at khub in the condensed 
phase for lo 7^ 00 in SCA model as in ZRP. However 
as shown in the main plot of Fig. 6 (a), m k increases 
linearly in k up to khub and there is no jump unlike in 
ZRP. In Fig. 6 (a), we plot m& versus k for p = 3.0 and 
lo = 1.0 which corresponds to the condensed phase of 
both RN and SFN with 7 = 4.3. We also measure mk in 
the fluid phase (p = 0.2, a; = 1.0) of the both networks 
(not shown), and confirm mk ~ k without any jumps. 
For SFN with 7 = 2.4, m k also shows the linearity in k 
for p — 3.0, lo = 1.0. Our simulation results imply that 
the relation ~ k is valid for any p > in RNs and 
SFNs with 7 > 2 for finite diffusion rate. The diffusion 
of masses indeed results in the linearity of nik in k for 
any lo 7^ 00. 

In the steady state, as masses can perform random 
walks with a finite rate, the mass of a hub node (vrihub) 
diffuses to different nodes with the probability of being 
at node i given as Eq. (5) [l^. The inset of Fig. 6 
(b) shows the snapshot of the mass distribution of nodes 
with degree k for a single sample at time t — 4 x 10 5 
on SFN with 7 = 2.4 for p = 3.0 and lo = 1.0. As 
shown, there is a peak which may be formed at khub as in 
ZRP. However the peak of mass m* is not always located 
at k} m b but diffuse around nodes according to P°° (Fig. 
6 (b)). Hence by taking average, the peak soaks into 
the average mass mk unlike in ZRP where all samples 
have the peak at khub- To see this more explicitly, we 



derive the relation, ~ k, based on the assumption 
that the diffusion (the random walks of masses) is the 
only relevant physical factor to decide P(m) in the steady 
state. 

First, we consider the average total mass Mk of nodes 
with degree fc defined as 



E 



mPo^ra, fc) 



(7) 



where Poo(m, k) is the probability of finding a random 
walker with mass m at degree k in fc-space in the steady 
state. As mass distribution P{m) in the steady state is 
independent of k, we have Poo(m, k) — P(m)PL x . P£° 
is the probability of finding a random walker at a node 
with degree k on the network. Then using P°° of Eq. JSJl, 
we write P£° as 



N 



i=l 



kNP(k) 



(8) 



where P(k) is a degree distribution. From Eq. (JJJ and 
©, we have 



M k 



kNP{k) 



mP(m) 



,m=0 



(9) 



So M fc = mkP(k)/ < k >~ kP(k). Then the average 
mass of a node with degree k is given as 



m k 



M k 
NP(k) 



k/N. 



(10) 



As expected, mk scales as mk ~ k for any degree distribu- 
tion P{k). As Eq. 1)10(1 is valid for any density p > 0, mk 
does not undergoes the condensation transitions unlike 
P(m) of Eq. (1). Therefore the steady state distribution 
of mk is determined by diffusion of masses rather than 
chipping of unit mass unlike in ZRP. 



VIII. SUMMARY 

We investigate the condensation phase transitions of 
symmetric conserved-mass aggregation (SCA) model on 
networks. In SCA model, masses diffuse with unite rate, 
and unit mass chips off from mass with rate lo. SCA 
model undergoes condensation phase transitions via the 
competition between diffusion and chipping processes 
0. 

First we consider the lo ^ 00 case. On random and 
scale-free networks of 7 > 3, SCA model undergoes the 
same type of condensation transitions from fluid phase 
into condensed phase as in one dimensional lattice of Rcf. 
13 . However unlike on regular lattices, the critical line 
on the networks depends on the network structures. On 
the other hand, on scale- free networks of 7 < 3 where one 
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or several nodes have a finite fraction of degrees, an infi- 
nite aggregation with exponentially decaying background 
mass distribution always takes place for any nonzero den- 
sity, so no phase transitions occur for oj =/= oo. The 
condensation and exponential mass distribution of small 
masses come from the generic inhomogeneity of network 
structure of SFNs with 7 < 3. However we are not able to 
numerically show the existence of the condensation phase 
at zero density limit for 7 < 3 due to the small world na- 
ture of networks |5j- Instead we numerically study the 
survival probability of a particle in pair annihilating ran- 
dom walks, so-called one lamb and one lion problem on 
networks. 

For the formation of an infinite aggregation of masses 
at zero density limit, unit mass chipped off from the in- 
finite aggregation should aggregate again with the ag- 
gregation within finite time interval. We numerically 
show that in thermodynamic limit, the survival prob- 
ability S(t) of a lamb (unit mass) is finite on random 
networks and scale-free networks (SFNs) with 7 > 3, but 
exponentially decays to zero with finite life time on SFNs 
with 7 < 3. Based on the finite life time of a lamb on 
SFNs with 7 < 3, we indirectly confirm the incomplete 
condensed phase for any p > on SFNs with 7 < 3. 

At to = 00, SCA model corresponds to zero-range pro- 
cess (ZRP) with constant chipping rate pj. ZRP with 



constant chipping rate on SFNs was shown to exhibit 
complete condensation [||. We also numerically show 
that the complete condensation takes place for any p > 
on random networks by measuring average mass of a node 
with degree k. Hence the complete condensation always 
takes place for any nonzero density in ZRP with constant 
chipping rate on random and scale- free networks. 

Finally we investigate the behavior of the average mass 
of a node with degree k, rrik in the fluid and the con- 
densed phase. In ZRP with constant chipping rate, rrik 
linearly increases with degree k, and jumps to the to- 
tal mass of the system at hub degree khub 0- Hence 
in SCA model, rrik is expected to show such a jump in 
the condensed phase for uj ^ 00. However rrik linearly 
increases up to khub without jumps in both the fluid and 
the condensed phase. We numerically confirm the linear- 
ity of rrik, an d also analytically show rrik ~ k with the 
assumption that the diffusion is the only relevant factor 
in the steady state. Therefore the steady state distribu- 
tion of rrik is determined by diffusion of masses rather 
than chipping of unit mass unlike in ZRP. 
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10148-0 from the Basic Research Program of KOSEF. We 
also thank Prof. Soon-Hyung Yook for critical reading of 
our manuscript. 
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FIG. 1: P(m,t) of m = 1 and 10 for p = 0.1 and w = 1 
on networks of A*' = 10 4 . Each panel corresponds to RN (a), 
SFN with 7 = 4.3 (b) and 7 = 2.4 (c). The vertical dotted 
line denotes the characteristic time r m . 
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FIG. 2: P(m) on RN (a) and SFN with 7 = 4.3 (b) at lo = 1.0. 
In each panel, the solid and the dashed line correspond to 
P(m) of p = 3.0 (condensed phase) and 0.2 (fluid phase) 
respectively. 




FIG. 3: The p-uj phase diagram of RN (a) and SFN with 
7 = 4.3 (b). In each panel, the solid line is the critical line 
of p c (u) = V 1 ^ + 1 — 1 on regular lattices of Ref. [l3j |. The 
dotted line between data points is a guide to the eye. The 
dashed line denotes ui = 00 line. 
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FIG. 4: P(m) on SFN with 7 = 3 (a) and 7 = 2.4 (b) 
at uj — 1.0. In each panel, the solid and the dashed line 
correspond to P(m) of p = 3.0 and 0.2 respectively, (c) the 
typical phase diagram of SFNs with 7 < 3. The dashed line 
denotes the line of ui = 00. 




FIG. 5: Semi-log plot of the survival probability S(t) on SFN 
with 7 = 4.3 (a) and 2.4 (b) of various size N up to 10 6 . (c) 
The average life time T on RN and SFN with 7 = 4.3. (d) T 
on SFNs with 7 < 3. The solid line between data points is a 
guide to the eye. 




FIG. 6: (a) The main plot shows mu of Eq. (10) for p — 3.0 
and u) = 1 on RN, SFN with 7 = 4.3 and 2.4 in condensed 
phase. The inset shows mt at u) = 00 on RN. (b) Mk of a 
single sample on SFN with 7 = 2.4 at different time for p — 3 
and w = 1 of the condensed phase. The main plot only shows 
the maximum mass m* of Mk of a single sample at different 
time steps, and the inset shows the snapshot of the M k at 
t = 4x 10 5 . 



